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The zitterbewegung being proportional to sin(et), it is depicted as the motion of electron from 
the positive energy state to that of the negative energy and vice versa in the neighbourhood of 
Dirac point. Since such transition involves crossing the Dirac point, the same may be the possible 
realisation and measure of the zitterbewegung. In this respect, we draw a similirity between the 
zitterbewegung and the Adler-Bell-Jackiw anomaly in (1+1) QED. In a time dependent perturbation 
theory, the embodiment of zitterbewegung is envisaged in the context of graphene. 

PACS numbers: 03.65.Pm, 72.80.Vp 



I. INTRODUCTION 

Electron in graphene, a monolayer of graphite with 
honeycomb lattice structure, exhibits relativistic nature 
as evident from the observed linear dispersion of energy 
and the zero energy state causing anomalous quantum 
Hall effect [3-|3|. While in the reciprocal lattice space 
there are two valleys of opposite momentum, K and K' , 
the space of honeycomb lattice consists of two inequiva- 
lent sublattices, A and B. In the low energy limit, the 
electron in graphene is described by the Dirac equation: 
the two valleys are analogue of energy or momentum and 
the sublattices that of spin. Because of the relativistic- 
like character of electron in graphene, the phenomenon of 
zitterbewegung of electron [5|-l7| has been shown renewed 
interest to unravel its the presence [§-fl4|. 

Schrodinger predicted the zitterbwegung of electron in 
the Dirac descrption by following the equation of motion 
for the position operator, r: 



dt 



(1) 



since [H, r] = —ihca where H = ca ■ p + /3mc 2 is free 
particle Dirac Hamiltonian, a and f5 the Dirac matrices, 
p the momentum, m the electron mass and, c the velocity 
of light. The relation in eq. ([T]) is independent of the 
electron mass. This is in contrast with the corresponding 
non-relativistic case 



dx 
~dt 



P_ 

m 



(2) 



with the Hamiltonian H = p 2 /2m. The eq. ([T]) is the 
well-known statement that in special relativity the ve- 
locity and momentum are fundamentally different. From 
the eq. ([T]), the position at time t is |5| 



= r(0) +c 2 H- 1 pt 



-ihc m H~ 1 e~ 2lHt ' h 



(3) 
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where t]q = a — cH~ 1 p which is obtained using the equa- 
tion of motion for a in the Heisenberg picture and r(0) 
the position at t = 0. In eq. ([3]), the second term is 
the average value while the third is instantaneous which 
is oscillating rapidly between the positive and negative 
energies. This oscillation is known as the zitterbewegung 
of electron [5[- The oscillatory behaviour is interpreted 
to be arising due to interference of positive and nega- 
tive energy states [Hj]. Let us note that in Dirac theory, 
an initial state of mean momentum and spin (+p, +j) 
given by a Gaussian wavepacket will get evolved in time 
t into a state that contains an admixture of (—p, ±j) and 
(+p, +j) if the wavepacket width is smaller than or equal 
to the Compton wavelength of the electron. Even if the 
required confinement is achieved, the oscillatory compo- 
nents is too small to be predictable. Therefore, one looks 
for a process being an embodiment of the zitterbewegung. 

This article attempts to depict the zitterbewegung phe- 
nomenon as the motion of the electron involving transi- 
tion between the positive and negative energy states near 
the Dirac point (zero energy level) in energy-momentum 
space and draw a similarity to the chiral or Adler-Bell- 
Jackiw (ABJ) anomaly [16| in (1+1) diemensional QED. 
Since the transition involves crossing the dirac point, the 
same may be a possible realisation and measure of zit- 
terbewegung. We formulate such embodiment of zitter- 
bewegung in a time dependent perturbation theory. The 
formulation is extended to graphene. The claim of simi- 
larity is limited to an anomaly-like realisation of the zit- 
terbewegung. It does not imply one can be given in terms 
of the other. 



II. ZITTERBEWEGUNG VS. ABJ ANOMALY 

Let us consider the motion of the electron in the neigh- 
bourhood of the Dirac point. An electron state initially of 
momentum and spin being (p, s) crosses the Dirac point 
(the zero energy level) from its positive neighbourhood 
into its negative neighbourhood where it acquires new 
momentum and spin (— p, r). Continuous occurence of 
such transition of one state of (p, s) to another of (—p, r) 
and vice versa in course of time results in the oscillation. 
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Each transition involves crossing of the Dirac point. This 
happens at a finite interval of time. This oscillation is the 
zitterbewegung of electron. The periodicity of the oscilla- 
tion gives the number of the Dirac point crossing. That 
is the measure of zitterbewegung. The transition does 
not necessarily change the spin but the momentum. 

The zitterbewegung part of the Dirac current as given 
in eq. (|A2[) (The notation follows the appendix.): 

(^7 M V) Z « sin(2e< + <p) (4) 

where the subscript Z means the zitterbewegung. The 
rhs in the equation above is the only time dependent 
part of the current. The period of oscillation depends on 
the range of values which e, <f> and t take. This means the 
there is a restructuring of e. This is the feature of zit- 
terbewegung, despite the fact that it is kinematic effect. 
It is almost impossible for the restructuring to happen of 
QED vacuum without any perturbation. For, a pertur- 
bative analysis is proposed in section IIIII 

This picture of zitterbewegung naively appears similar 
to the chiral anomaly in (1+1) QED which is a quantum 
effect though. 

Let us recall the idea of the ABJ anomaly in (1 + 
1) dimensional QED [l7|. The mechanism of the ABJ 
anomaly relates to the production of Weyl fermions in 
the presence of an external field, A\ . The so-called Weyl 
fermions are massless spin i particles described by a par- 
ity non-invariant Hamiltonian. They are of two compo- 
nents, left-handed (RH) and right-handed (RH). In the 
absence of At, both the LH and RH fermions meet at 
the point of degeneracy. On slowly applying the exter- 
nal field, A\ , the degeneracy is lifted and the LH fermion 
moves down below the degenerate level while the RH 
fermion moves up in the energy-momentum space. For 
a particular value of A\ both LH and RH fermions meet 
again at the point of degeneracy. On further application 
of A±, the degeneracy is lifted again and the LH fermion 
moves up while the RH fermion moves down. It should 
be understood as restructuring of energy levels or change 
of LH fermion into RH and vice versa while crossing the 
point of degeneracy as shown in Figure [T] The result is 
the nonconservation of axial charges. That is the ABJ 
anomaly. The nonconservation of axial charge is seen 
as the movement electrons from the neighbourhood of 
LH degeneracy point to that of RH in energy-momentum 
space, without net production of electrons, in the pres- 
ence of external field [18| • 

Although the two phenomena are different, the naive 
similarity arises at the level of restructuring of the energy 
levels which corresponds to the transition of a state from 
the positive energy to the negative. It is interesting to 
note the following: the chiral anomaly (a) occurs due to 
mass term in the absence of external field (Fujukawa der- 
vation [l9| ) , (b) chiral current is not conserved and (c) is 
sensitive to dimensionality, whereas the zitterbewegung 

(d) involves (although implicitly) spin-orbit interaction, 

(e) implies non-conservation of spin current and (f ) is in- 
dependent of dimensionality. However, effective dimen- 



E 




FIG. 1: Chiral anomaly: shifting of energy levels or transition 
LH (solid line) and RH (broken line) fermions as a function 
of Ai . Spatial length is set to be unity. 

sional reduction occurs when the system is under external 
magnetic field. While zitterbewegung is a kinematic ef- 
fect as a consequence of incorporating special relativity 
in quantum mechanics, the chiral anomaly is a dynami- 
cal phenomenon as a consequence of second quantisation. 
From the point of their similarity, we may get an insight 
of zitterbewegung of electron. 

III. ZITTERBEWEGUNG CURRENT 

In order to show how an initial state of momentum p 
and spin s can make a transition into a state of momen- 
tum — p and spin r in presence of slowly varying exter- 
nal potential, we consider a time dependent perturbation 
theory. Embodiment of the zitterbewegung is elucidated 
in the resulting S matrix. 

For a free massless electron in presence of a slowly 
varying external potential, V^x, t). 

crfp^{x) = -erfV^x, t)*(x) (5) 

The solution to the above equation is, in mode expansion, 

H>(x) = ip(x) 

+ J2 [<s(tKe^ + a- r (t)v r p e-^] (6) 

p,s,r 

where if)(x) is the solution of free particle Hamiltonian 
and can be taken as well defined initial state of (p, s) and 

a+ s (t) = -e| J d 3 x' J dt'u^V^x'^ixy-^) 

a; <s (t) = ei Jd s x' J + dt'v s p ^V^x',t')i;(x')e ipx '(8) 

where E is the eigen energy. Of the terms in the second 
line of eq. ([6]), the first term is the probability of the 
initial positive state remaining the same with ti = — oo 
while the second one is the probability of a positive en- 
ergy state making transition into a negative energy state 
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with momentum —p and spin r in presence of external 
field. To make this statement hold, the t' integration in 
eq. (jHJ has to be taken between finite time interval as 
indicated. Taking one of the limits to be infinity implies 
creation or annihilation process. The equations in ([7]) 
and (jHJ are S'-matrices of first order in e expansion in 
presence of external potential. 

With the limits of integration specified in equation © , 
the t' integration turns out to be: 



_ iEt , = 2sin(£tt) 
E 



(9) 



The scattering matrix corresponding to eq. © is 



E . ' . (10) 
The S-matrix in the above equation can be shown to 

be containing the current driving the transition of the 

positive energy state to the negative energy one. For, 

we take derivative of S u -w with respect to the external 

potential V^(x, t). 

Since we have chosen the potential to be applied adia- 

batically, we have 



V^x,t) = e- A %(a:) 



(11) 



where A is chosen to be a small parameter signifying 
switch on and off the potential. The variation of V(x, t) 
with respect to time is negligibly small. Therefore, 
dV = V df(x) if V(x) = Vof(x) and otherwise dV = dV . 
The current is obtained as 



E 



dS u 



(12) 



A simple comparison of the current term given by eq. 
(|12|) with the one in eq. (jlj leads to identification of 
current given by the scattering matrix, eq. (|12[) . is the 
zitterbewegung current. 



IV. GRAPHENE CASE 

We would like to extend the discussion in the previous 
section to graphene. 

The Dirac equation for a massless electron in graphene 
in the presence of external magnetic field is 



= 



(13) 



where Vf is the Fermi velocity of electron in graphene 
which is 10~ 2 c, the Greek indices p = 0,1,2, = 
(A t ,A x ,A y ) — (0, 0, Bx) and p x is expressed in units of 



H/Ib where Is 



he 



The energy spectrum is 



eB 



E n = ±hv f \j—{2n + 1-0-3) 



(14) 



The n = energy levels are: 



En 




(15) 



The above defines the Dirac point in the presence of ap- 
plied magnetic field. The solution of eq. ([TBI is 



ipn(x) = exp[ip y y] exp 



eB (x — x ) 2 
lie 2 



(16) 



where xn = j§k y the extension of the wavepacket, and 
w 3 , j = 1,2,3,4 are normalised spinors. 

In the presence of magnetic field the electron states 
are Landau quantised. The electron is confined into a 
"strip" of width £0 . Now, We look at the motion of the 
electron in presence of external electric field, Vq applied 
adiabatically, parallel to the magnetic field. Then, the S- 
matrix obtained from eq. © is, after the t' integration, 

S = § / d?x'vL p l°V i>n{x')e-^ s ' sm(Et) (17) 



From eq. (|17|) following eq. (|T2j is the charge den- 
sity, j°(t). Obtaining the eigen energy E = E(B,Vo) 
for n = 0, we can find the oscillation and dermine the 
energy restructuring in course of time. As mentioned be- 
fore in section \H\ the measure of zitterbewegung can be 
obtained. 



V. CONCLUSION 

Although zitterbewegung is a kinematic phenomenon, 
it is an spin-orbit coupling in an unusual way. The spin- 
orbit interaction involves both positive and negative en- 
ergy states. It can be turned into a dynamical when an 
electron is subjected to an external field. Its signature 
is the crossing of Dirac point. At this level, its simi- 
larity with ABJ anomaly arises as discussed. Despite 
the similarity, they are not related. However, given the 
understanding of the chiral anomaly over decades, the 
anology may be useful to have better understanding of 
zitterbwegung. In graphene, the spin-orbit interaction 
is the pseudospin-orbit interaction. Pseudospin is given 
by the sublattice states. In the desription of graphene 
with chiral symmetry, the chiral degree of freedom is re- 
lated to the sublattice states. The chiral symmetry and 
its breaking in graphene have become interesting aspects 
of graphene recently [2314261 ] . It is interesting to note 
the possibility of chiral anomaly induced charge chang- 
ing transition in graphene (27j . In view of the similarity 
between the zitterbewegung and the chiral anomaly, we 
may look for a possible dynamical description of zitter- 
bewegung and realisation in graphene. 
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Appendix A: Appendix 

Given the wavefunction for a free massless Dirac par- 
ticle 

= E E J^ c P^ s (p)e ipx + 

p s=l,2 

E E \/Y s d irV r (P>- iPX (Al) 
p r— 3,4 
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